Dynamical Mean-Field Theory of Electron-Phonon Interactions in Correlated 
Systems: Application to Isotope Effects on Electronic Properties 
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We use a recently developed formalism (combining an adiabatic expansion and dynamical mean- 
field theory) to obtain expressions for isotope effects on electronic properties in correlated systems. 
As an example we calculate the isotope effect on electron effective mass for the Holstein model as 
a function of electron-phonon interaction strength and doping. Our systematic expansion generates 
diagrams neglected in previous studies, which turn out to give the dominant contributions. The 
isotope effect is small unless the system is near a lattice instability. We compare this to experiment. 
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I. INTRODUCTION 



In solid state physics the term "isotope effect" has 
come to mean a dependence of an electronic property 
on an ionic mass. The isotope effect on the supercon- 
ducting transition temperature Tc of conventional super- 
conductors was important evidence fap, the role played 
by phonons in the pairing mechanism. EJn In conventional 
metals the isotope effect on electronic properties other 
than the superconducting Tc is negligible. The Migdal- 
Eliashberg (ME) theory of electron-phonon coupling in 
metalsQ explains this as follows: most electronic proper- 
ties are determined by processes occurring on the scale of 
the electron kinetic energy t. Interaction with phonons 
affects electronic properties only on the scale of a typical 
phonon frequency ujq, which is much less than t. Thus 
isotope effects are generically expected to be small, of 
the order of the adiabatic parameter j — ujo/t <^ 1. ME 
theory essentially retains terms only of order 7°; there- 
fore, most isotope effects are beyond the scope of this 
theory. The only exception is the superconducting tran- 
sition temperature Tc = u;oexp(— 1/[A — ^*{ujo)]) where 
the phonon frequency enters as the upper cutoff of the 
logarithmic divergence in the pairing interaction and the 
lower cutoff of the logarithmic divergence in the Coulomb 
pseudopotential n* . 

Recent observations in several classes of "strongly cor- 
related" materials (high-temperature superconductors, 
"colossal magnetoresistance" manganites, and alkali- 
metal doped Ceo) of large isotope effects on eh 
tronic properties, including electron— efective masse 
and supercgnducting,B a magnetic, t^TEj and charge 
orderingliim transition temperatures pose a fundamen- 
tal challenge to this understanding ancL gall -for a the- 
ory that goes beyond ME. Some authorsErlZlll3'E£l have at- 
tempted to relate their experimental data to formulas 
derived for the case of a "polaron": a single electron 
interacting with lattice deformation. All experimentally 
relevant systems, however, have a metallic density of elec- 
trons, of order one per unit cell, so thei applicability of 
"polaron" formulas is not clear. Otherslla have considered 



the first corrections to ME theory for the electron mass 
using standard diagrammatic methods but have not con- 
sidered the feedback effects of electrons on phonons and 
have presented results that cannot easily be extended to 
correlated systems. _ 
Very recently, we have introduced a new metiiod,Ej 
combining the dynamical mean- field (DMF) theorjCa and 
an adiabatic (small-7) expansion, for studying electron- 
phonon interactions in systems with arbitrary electronic 
correlations. Here we use this method to calculate 
the electron effective mass m* due to interactions with 
phonons. Our analytic result includes order-7 quantum 
lattice fluctuations that were neglected by ME. These 
fluctuations give rise to a non-zero isotope effect on to*, 
which is small unless the system is sufRciently close to 
a polaronic instability. Within our model, the isotope 
effect is negative near half filling and positive away from 
half filling, which is due to the competition of lattice and 
density fiuctuations. 

The paper is organized as follows. In Sec. |l^ we in- 
troduce the electron-phonon lattice Hamiltonian and its 
associated effective phonon action within the local DMF 
formalism and recall the main ideas underlying the adia- 
batic expansion. In Sec. Ill we calculate physical quanti- 

elec- 



In Sec 

ties "beyond ME" : the Luttinger-Ward functional 



tron self-energy, and isotope effect on electron effective 
mass. In Sec. IV we compare our results to recent ex- 



perimental findings and draw some general conclusions. 

II. MODEL AND METHOD 

Wc study a general tight-binding based Hamiltonian 
H = Hci + i?ph + Hei-ph, where 

Hc\ = - ^ U^j {cl^Cja + c]c,C.i„) - N^JLTl + Hcc, (1) 

Hp,. = \Y.^Mxj + Kx% (2) 

i 

Hci-ph ^ g^Xi{n^ - n). (3) 
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The operator cj^ creates an electron with spin a on lattice 
site i. The mean density n = (l/N) claCicr (where N 
is the number of lattice sites) is fixed by adjusting the 
chemical potential /z. Electron-electron interactions Hcc 
are not explicitly written. The operator Xi measures the 
ionic displacement at site i. From spring constant K and 
ion mass M we can define the characteristic oscillator fre- 
quency Wo = {K/My^'^, which we take to be dispersion- 
less (Einstein model). The electron-phonon interaction 
couples the phonon displacement Xi to the electron den- 
sity n-i = J2a '-L'-ic the same site. We define (xi) =0 
to be the equilibrium phonon displacement for a uniform 
electron distribution. 

In DMF theory the properties of TJ—tpay be obtained 
from the solution of an impurity modelEj specified by the 
action S[c, c, x, a] — So[x] + Scc[c, c, a] + Si[c, c, x, a], with 



(4) 



Si[c, X, a] — — y ^ Cna-Cna-an +.9 ^ ] CngCn+k.a-Xk ■ (5) 
n a nka 

The impurity electron is represented by Grassmann fields 
Cno- , Cna , which depend on odd Matsubara frequencies 
u)n = (2n + 1)'kT and spin a =1, i- The impurity phonon 
is represented by a bosonic field Xk, which depends 
on even Matsubara frequencies Uk = 2fc7rT. Electron- 
electron interactions (arising from -ffoe) are described by 
5cc- The local electron Green function is defined by 
Gioc[fl]n ~ (5 In Z[a]/5a„ where 



Z[a] 



[dcdcdx] exp —S[c, 



c, X, a\ 



(6) 



is the partition function. The mean field function a 
(which contains information about the nonlocal physics) 
is fixed by the self-consistency condition 



Gioc[a]Ti — / dek 



fi - E[a]„ - ek 



(7) 



which equates Gioc [a] and the momentum integrated lat- 
tice Green function. The momentum integral has been 
converted to an energy integral using the lattice density 
of states p(efe). 

We integrate out the electron fields and work with 
an effective phonon action S'[a::, a] = Sq[x] + 
S'i[2:,a] — — In J[c?c(ic] exp(—S'[c, c, x, a]). In practice S 
may depend on additional auxiliary fields (such as spin 
and charge fluctuation fields), which have to be aver- 
aged over. We then use the crucial facia that the scale 
cjQ = {K/My^'^ on which the phonon fields x vary is 
much smaller than the scale t (bandwidth or interaction 
scale) on which electronic quantities (such as G and a) 
vary so that an expansion is possible in the "adiabatic 
parameter" j — LUo/t <^ 1. In a first step we formally 
expand Si in powers of x about a; = (corresponding to 
a conventional metallic state with no lattice distortions), 



ig0 a;" 



v"/2 



In a second step [to be performed 
below, in the analysis leading up to Eqs. (p^)-(p^)] the 
phonon vertices of the action may be evaluated via a low- 
frequency expansion. 



III. PHYSICAL QUANTITIES 

In the following we will not directly work with the ef- 
fective phonon action but with the Luttinger-Ward func- 
tional derived from it. In general (j) = (jf° + (fp^ is 
defined as the sum of all vacuum-to-vacuum skeleton di- 
agrams. Within the DMF appxpximation it is related to 
the local free energy fiimp viaE3 



m 



(lnG„ 



(8) 



The electron self-energy E = -I- is a saddle point 
of r^imp. This implies E ~ 6(j)/5G. Since the ther- 
modynamic potential r^imp [given by Eq. (||)] and the 
partition function Z [given by Eq. (H)] are related via 
exp(— /Siliinp) = Z we may obtain (fp^ from an adiabatic 
expansion of Z . To order 7^ we find 



d^^''[G] = W\nDl^ 



^ ^i[G]k,k' -k-k'DkDk' 



kk' 



Y.h[G]lk',k-k'DkDk'Dk-k', 



(9) 



kk' 



where Dk = {1 + {ujk/ojof - Af2[G']fc,_fe) ^ is the 
phonon Green function and '^n[G]ki,...,kn dimension- 
less phonon vertices whose explicit form depends on Sec- 
The phonon fields gx ^ x were rescaled to have units 
of energy. The parameter A = g"^ /{Kt) is the electron- 
phonon interaction strength. The Feynman diagrams 
corresponding to Eq. (|^) are shown in Fig. |^. 

We use Eq. (^ to formally calculate S to 0(7^). 
Following the "adiabatic expansion" outlined above we 
may further simplify the resulting expressions by evalu- 
ating the r„ (whose dominant contributions come from 
frequencies of order t) via a low-frequency expansion. 
To the given order in 7 we may approximate the cu- 
bic and quartic vertices by their static value (which we 
write Fs and r4) and the phonon Green function by 
Dk « (1 - A/Ac + (wfeM)' + A7ap|wfe/wo| + 0(7'))"', 
where Ac — ^2[G]q \ is a critical interaction strength and 

Qfp is a damping parameter. The static r2 vertex acts 
as a phonon self-energy and renormalizes the expansion 
parameters to 7 = 7(1 — A/Ap)^/^ and A = A/(l — A/Ac). 
If A — > Ac from below then A — > 00. This "polaronic in- 
stability" was analyzed in great detail in Ref. ^ In the 
following we assume A < Ac. The self-energy to 0(7^) 
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can be written as T,^^ = Sia + Sib + Sic + '^2a + S26 (see 
Fig. where 



Slaji — ^tT Dk Gn+k , 



(10) 



Sib„ — A'^f4T^ ^ D^iDfe/G„+fe, (11) 

kk' 

Sicn = 2A^f3T^ ^ _D^._Dfe/Dfe-fe'G„+fc, (12) 

kk' 

S2ari = ^t^T^'^^DkDk'Gn+kGn+k'Gn+k+k' , (13) 
fcfc' 

S26n = 2}?tt''iT'^'^^DkDk'Dk-k'Gn+kGn+k' ■ (14) 

fcfc' 

The 0(7) one-loop diagram Ei^ is at the basis of ME 
theory. The 0(7^) two-loop diagrams Sif,, Sic, S2a, and 
S2b represent lowest-order corrections to ME. Apart from 
Sic they all arise from vertex corrections of the electron- 
phonon coupling. The diagrams Sib, Sic, and S26 were 
overlooked in Refs. [iq and 
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We now turn to the electron mass enhancement from 
interactions with phonons, defined by m*/m |ph = 1 — 
5SPh(w)/a(ii/)|^=o ^l + \la + -.- + Mb, with SPh(izy) 
written at T = 0. We will calculate m*/m to 0(7). 
Since the frequency derivative is of 0{j~^) we need S 
to C(7^) as above. We first calculate the mass enhance- 
ment Ala from the basic one-loop diagram Sia(ij/) = 
Xt J duj/ {2'K)G{iv + iu])/{l - A/Ac + (w/t^o)^)- In gen- 
eral Iml]iafl(z/) is of order A7 and is zero for \v\ < wq. 
From this and Eq. (@) it follows that G(w) = G'"=(ii/) -I- 
G'=°'>(w) where 

= + A^)' = -*^sign(j.)p(A.) 



2t 



(15) 



at low frequencies, for Sec — 0, and neglecting terms of 
0(7). Weak electron-electron interactions renormalize 
the scale t in Eqs. (|l^), but their general form remains 
the same. We thus find 



AiQ = \tp{fi) - -A7, 



(16) 



where the leading term comes from the derivative acting 
on G*^"^, producing a delta function 6{uj), and the 0(7) 
term comes from G'"'^. The phonon damping ~Jw| in D 
does not give rise to any 0(7) corrections to (16f). The 
two-loop diagrams can be analyzed in a similar fashion. 
We just give the results: 

Aifc = ip7tp(M)f 4, (17) 

Ale = h^^tp{f,)fl (18) 

A2a = -I'X^ltpip) {i-Kh-'pipf - 5{p/2tf) , (19) 
A2fc = ~X^^tp{p){p/2t)r3. (20) 



In principle our adiabatic expansion is based on the 
smallness of the parameter 7 alone, with no assumption 
being made about A. However, 0(7") terms in the ex- 
pansion are of the form A^t"" ~ (1 — A/Ac)"/^~™ with 
m > n and become arbitrarily large if A Ac. We there- 
fore need to set the range of parameters A and 7 where 



the expansion is "\ 
First we notice 



id. We propose two different criteria, 
that at half filling m*/m as a func- 
tion of A has a local maximum at Amax 1 and then goes 
to — (X) as A — > cxD. This unphysical behavior comes from 
the fact that T3 vanishes at ^ = and that the dom- 
inant term Ak, ^ A'^7r4 to this order is negative. For 
< \p\ < pi ^ 1, i-e. sufficiently close to half filling, the 
local maximum at Amaxi can still be defined, but now 
m* I'm — > -\-oo as A — > 00 because the dominant term 
Ale ~ A'*7r| is positive. The upper bound pi depends 
weakly on 7, as shown in the first two columns of Table 
^ We thus suggest that for given 7_and \p\ < pi{j) our 
adiabatic expansion is valid if A < Amax 1 • In Table | we 
give Amax 1 and the corresponding unrenormalized quan- 
tity Amaxi at half filling {p = 0), for a few representative 
values of 7. 

Farther away from half filling, i.e. for \p\ > pi, the 
local maximum in m* /m ceases to exist and we need 
a different criterion Amax 2, which we define as follows: 
for given 7 we need A < Amax 2 small enough so that 
the ratio of the 0(7) terms and the 0(1) term in m* jm 
be less than p < 1. We fix the parameter p by setting 
Amaxi = Amax 2 at /J, = 0: the result is p w 0.25, which is 
compatible with the 0(7) terms being a small perturba- 
tion. In Table |l| we give Amax 2 and Amax 2 at half filling 
and at p/{2t) = 0.4, for some values of 7. 

We now return to the effective mass. In the fol- 
lowing we assume Sec = and a semicircular den- 
sity of states p(efe) = l/(27rt2)(4i2 _ £2^1/261(4^2 _ ^2)^ 

for which the vertices can be evaluated exactly: r2 = 
8[1 - (Ai/2i)2]3/2/(3^), fa = ^{pl2t)[l - {p/2tfY/^/^ 

and f4 = -16[1 - {p/2tfY/^[l - 6(/i/2t)2]/(157r). In 
Fig. ^ we plot m* jm, as a function of A for three differ- 
ent values of 7. At half filling (/i = 0, main graph) we 
plot up to Amax 1 (cf . Table ^ . Far away from half filling 
/(2<) = 0.4, inset) we plot up to Amax2(0.4) (cf. Table 



The isotope effect on electron effective mass is defined 
by Qfm* = — c?lnm*/(iln Af . Using M 7"^ it is conve- 
nient to rewrite this as am' — ^ / {2m*)dm* / d'y , so that 



Oi-tn* 



-iA7 + Aib + ... + A 



2b 



2(1 + \tp{p)) 



Oh'). 



(21) 



In the main graph of Fig. |j we plot am* as a function 
of A < Amax 2(0) (cf. Table at half fflhng and for 
different values of 7. The isotope effect is negative and 
its absolute value is limited by 0.125 = p/2. In the inset 
of Fig. H we plot am' at p/{2t) — 0.4, for the same val- 
ues of 7. For 1.5 ^ A < Amax2(0.4), the isotope effect is 
positive. We explain this in the next section. 
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IV. DISCUSSION AND CONCLUSION 

We start with a brief analysis of S. The correct 
two-loop self-energy consists of four diagrams, three of 
which, the diagrams containing the vjectices T3 and 
were neglected by previous authors£3'c3 Fig. || and 
Eqs. (pi]), ( p^ suggest that the diagrams Sib and 
Sic can be absorbed in a 7-dependent Ac as follows: 
Aji(7) =t2 + {Ai Xfi + Bi f §) 7 + 0(72). If we as- 
sume T ^ ujq (which is appropriate for calculating m*) 
then Ai and Bi (which arise from the summation over 
Dk) are positive constants of order one. As pointed out 
in the previous section, the 0(7) term in Ac (7) changes 
sign when the system is doped away from half filling. The 
same is true for the mass enhancement A2a + A2b from the 
remaining two diagrams. We explain this as follows. Very 
close to half filling, low-frequency lattice fluctuations re- 
duce the average local ion displacement x. As a result, 
the electron-phonon coupling in Eq. (^) is less efficient, 
m* /m decreases relative to the value 1 + Ai^ correspond- 
ing to a completely static lattice, and the isotope effect 
is negative. This is the physics behind the terms Ai;, and 
A2a and the second part of Aia- It is reflected in the main 
graphs of Figs. |^ and ^. Away from half filling, on the 
other hand, the physics is dominated by density fluctua- 
tions, which have the opposite effect. They enhance the 
density-coupled interaction (^, increase m*/m, and lead 
to a positive am* ■ This is represented by the terms Aic 
and A2f, and in the insets of Figs. ^ and ^. 

We now comment on the subject of expansion param- 
eters, which has been the source of some controversy in 
the literature. The original ME articlestl and more recent 
worked indicate that A7 is the proper expansion parame- 
ter of a nonadiabatic theory, whereas other authors (see, 
e.g., Ref. p5|and references therein and, for a particularly 
clear discussion, Ref. ^2|) argue that the expansion breaks 
down if A exceeds a critical value of order one, irrespec- 
tive of 7. Within our method, it is possible to reconcile 
these two points of view. By working with renormalized 
parameters A and 7 we implicitly take into account the 
presence of the polaronic instability at Ac- The condi- 
tion A < Ac of Refs. ^ and thus appears as a natural 
limitation of our theory. If A is sufficiently close to the 
instability at Ac (which is possible even for A < Aniax(7), 
as can be seen from Tables || and ||) then the dominant 
mass enhancement terms Aifc ^ (1 — A/Ac)^^^^ (close to 
half filling) and Aic ~ (1 — A/Ac)"^/^ (away from half 
filling) act by shifting Ac to higher and lower values, re- 
spectively, as shown in the previous paragraph. This 
argument can immediately be extended to higher orders 
in 7, noting that the odd vertices r2n+i all vanish at half 
filling due to partcle-hole symmetry. We thus conclude 
that the effective expansion parameters of our theory are 
A7 close to half filling and ^^7 away from half filling. 
This generalizes the results of Refs. || and |2^. 

Next we compare our calculations to experiments. Ref. 
^ presented measurements of the oxygen isotope depen- 



dence of the Meissner fraction of LSCO high-temperature 
superconductors. The assumption that the dependence 
arose from an isotope effect on the carrier mass m** 
(and hence the penetration depth) implied « —0.5. 
Within our theory, such large values of am*" imply that 
the material is at or beyond the polaronic instability; 
however, this assumption may not be consistent with the 
observed reasonably good conductance of high- Tc mate- 
rials. On the other hand, Refs. ^ and ^ have demon- 
strated that in the "colossal magnetoresistance" mate- 
rial Lai_a;Caa;Mn03 (0.2 < x < 0.5) the isotope effect 
on low-temperature properties including the carrier mass 
(determined via specific heat) is very small, even though 
electron-lattice interactions are believed to be strong in 
this material and even though the isotope effect-on— the 
ferromagnetic transition temperature Tc is large.ll£rll3 In 
the manganites it seems likely that the very large isotope 
effects on Tc are associated with a_phase transition that 
is now believed to be first order;L3 we speculate that a 
similar phenomenon might explain the data on the high- 
Tc material. 

In summary, we have presented a detailed theory of 
isotope effects in models with strong electron-lattice cou- 
pling. We have determined the correct expansion param- 
eters and have found diagrams overlooked in previous 
works. Our results suggest that a large isotope effect 
on electronic properties is very difficult to obtain in a 
metallic systepj..,— Unlike previous treatments of nonadi- 
abatic effectsJl3'Eil our formalism can easily be general- 
ized to include electron correlations, which enter in two 
places: (i) vertex and self-energy corrections renormalize 
the phonon vertices r„ and (ii) irreducible particle-hole 
scattering vertices renormalize the electron-phonon cou- 
pling g. Both effects are believed to strongly suppress 
the electron-phonon interaction. An extensive analysis 
of problem (i) for r2 in the presence of local electron- 
electron interactions can be found in Ref. Progress 
on (ii) will be presented in a future publication. 

We thank S. Blawid and R. L. Greene for useful dis- 
cussions. We acknowledge NSF DMR00081075 and the 
University of Maryland-Rut gers MRSEC for support. 
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FIG. 1. Luttinger-Ward functional to 0(7^). A wavy line 
stands for the phonon Green function D. The shaded circles 
stand for the phonon vertices T2, Fa, and r4. 
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FIG. 2. Self-energy diagrams to 0(7^). A straight line 
stands for G, a wavy line stands for D. The phonon vertices 
Fa and r4 are represented by shaded circles. 




FIG. 3. Electron mass enhancement m* /m due to inter- 
actions with phonons, as a function of renormalized elec- 
tron-phonon coupling A and for different values of adiabatic 
parameter 7 = uJo/t. Main graph: half filling (^ = 0) and 
A < Amaxi(7)- Inset: iJ./{2t) = 0.4 and A < Amax2(7)- 
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FIG. 4. Effective mass isotope exponent Om* = — dlnm*/ 
d in M as a function of renormalized electron-phonon coupling 
A for different values of adiabatic parameter 7 = ujo/t, at 
T = and within the range of validity of the adiabatic ex- 
pansion. Main graph: half filling (/i = 0). Inset: away from 
half filling (M/(2t) = 0.4). 



TABLE I. Limiting values of A and A as a function of 7 
at half filling, using the first criterion (valid for < /ii) 
described in the text. 



7 



fJ.1 



An 



l(0) 



An 



i(0) 



0.05 

0.1 

0.15 

0.2 

0.25 



0.062 
0.079 
0.091 
0.100 
0.109 



10.8407 

6.2804 

4.4458 

3.4168 

2.7456 



1.0626 
0.9920 
0.9313 
0.8760 
0.8244 



TABLE II. Limiting values of A and A as a function of 7 
at half filling and at n/{2t) = 0.4, using the second criterion 
described in the text. 



7 


Amax 2 (0) 


Amax 2 (0) 


Amax 2 (0.4) 


Amax 2 ' 


0.05 


8.5620 


1.0356 


4.7589 


1.1579 


0.1 


5.2225 


0.9613 


3.7243 


1.0846 


0.15 


3.8793 


0.9037 


3.2578 


1.0412 


0.2 


3.1268 


0.8557 


2.9782 


1.0109 


0.25 


2.6369 


0.8143 


2.7874 


0.9879 



(0.4) 
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